*****************************************************************
********* Replication Files for Complementary Mediation *********
**************** Author: Elizabeth J. Menninga ******************
******** Accepted: July 2020; International Interactions ********
************ Stata Version 14.2 Used for All Analyses ***********
*****************************************************************

use "Complementary Replication Data.dta" , clear
set more off

************************************************
********* Codebook for Categorial Vars *********
************************************************

**Creating one categorical variable for interaction between cont_lev & stay_any
*gen cat_lev=.
*replace cat_lev=0 if cont_know_ord==0 & stay_any==0
*replace cat_lev=1 if cont_know_ord==1 & stay_any==0
*replace cat_lev=2 if cont_know_ord==2 & stay_any==0
*replace cat_lev=3 if cont_know_ord==0 & stay_any==1
*replace cat_lev=4 if cont_know_ord==1 & stay_any==1
*replace cat_lev=5 if cont_know_ord==2 & stay_any==1

**Codebook Note:
*cat_lev=0--> no soft lev
*cat_lev=1--> relational lev with 1 side & no staying power
*cat_lev=2--> relational lev with both sides & no staying power
*cat_lev=3--> no relational lev but staying power
*cat_lev=4--> relational lev with 1 side and staying power 
*cat_lev=5--> relational lev with both sides and staying power

**Codebook Note:
*coord_single=0--> single party
*coord_single=1--> uncoordinated MP
*coord_single=2--> coordinated MP

*************************************************
********* Descriptives Reported in Text *********
*************************************************

*medtype==2 NGO or individual
*medtype==1 IGO
*medtype==0 state

tab medtype1 
tab medtype2
tab medtype3 
tab medtype4 
tab medtype5 
tab medtype6 
tab medtype7 
tab medtype8


tab medcode1 
tab medcode2
tab medcode3 
tab medcode4 
tab medcode5 
tab medcode6 
tab medcode7 
tab medcode8

tab outcome3

sum gdp_ratio_thousands, detail

corr gdp_ratio_thousands cinc_ratio
tab stay_any
tab cont_know_ord
tab numofmed
tab coord_single


***************************************
********* Table 1 Regressions *********
***************************************
set more off
ologit outcome3 c.gdp_ratio_thousands i.cat_lev numofmed i.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)
estimates store m1
ologit outcome3 c.gdp_ratio_thousands i.cat_lev numofmed i.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5 currentpk, cl(conflictid)
estimates store m2

*Changing baseline for coord_single
ologit outcome3 c.gdp_ratio_thousands i.cat_lev numofmed ib1.coord_single fullconf_duration_10years currentpk intensity_level terr_conflict polity2_v5, cl(conflictid)
ologit outcome3 c.gdp_ratio_thousands i.cat_lev numofmed ib1.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)

*With Interaction
ologit outcome3 c.gdp_ratio_thousands##i.cat_lev numofmed i.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)
estimates store m3
ologit outcome3 c.gdp_ratio_thousands##i.cat_lev numofmed i.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5 currentpk, cl(conflictid)
estimates store m4
*Changing baseline for coord_single
ologit outcome3 c.gdp_ratio_thousands##i.cat_lev numofmed ib1.coord_single fullconf_duration_10years currentpk intensity_level terr_conflict polity2_v5, cl(conflictid)
ologit outcome3 c.gdp_ratio_thousands##i.cat_lev numofmed ib1.coord_single fullconf_duration_10years  intensity_level terr_conflict polity2_v5, cl(conflictid)

esttab m1 m2 m3 m4 using ologit_gdp_ratio.tex, b(3) compress ///
		se label replace cells(b(fmt(3) star)  ///
		se(par)) starlevels(* .05)  ///
		varwidth(30) collabels(, none) noomitted
		
		
***************************************
************* Model 1 & 3 Plots *******
***************************************
ologit outcome3 c.gdp_ratio_thousands i.cat_lev numofmed i.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)

********* Figure 2 ********

margins cat_lev, vsquish predict(outcome(0)) saving(out0, replace)
margins cat_lev, vsquish predict(outcome(2)) saving(out2, replace)

clear
use out0
twoway (bar _margin _m1 if (_m1==0)|(_m1==1)|(_m1==2), sort color(gs11) barw(.5)) ///
    (bar _margin _m1 if (_m1==3)|(_m1==4)|(_m1==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub _m1, sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel(0 "No Rel. Lev." 1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(0(.15).90, labsize(large)) ///
	xtitle(" ") ///
	ytitle("Pr(Failure)" " ", size(large)) ///
	graphregion(color(white)) ///
	legend(order(1 2) label(1 "Without Staying Power") label(2 "With Staying Power") cols(2) region(lcolor(white)))  ///
	saving(twoway_outcome0, replace)

clear
use out2
twoway (bar _margin _m1 if (_m1==0)|(_m1==1)|(_m1==2), sort color(gs11) barw(.5)) ///
    (bar _margin _m1 if (_m1==3)|(_m1==4)|(_m1==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub _m1, sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel(0 "No Rel. Lev." 1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(0(.15).90, labsize(large)) ///
	xtitle(" ") ///
	ytitle("Pr(Peace Agreement)" " ", size(large)) ///
	graphregion(color(white)) ///
	legend(order(1 2) label(1 "Without Staying Power") label(2 "With Staying Power") cols(2) region(lcolor(white)))  ///
	saving(twoway_outcome2, replace)

grc1leg twoway_outcome0.gph twoway_outcome2.gph, saving(cat_bar.gph, replace) ///
	note("Note: 95% Confidence Intervals reported.", size(medium)) ///
	graphregion(color(white)) pos(12)

********* Figure 3 ********
use "Complementary Replication Data.dta" 

ologit outcome3 c.gdp_ratio_thousands##i.cat_lev numofmed i.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)
margins, at(cat_lev=(0(1)5) gdp_ratio_thousands=(0 .15 .27)) vsquish predict(outcome(2)) saving(m3_out2_mean, replace)
margins, at(cat_lev=(0(1)5) gdp_ratio_thousands=(0 .15 .27)) vsquish predict(outcome(0)) saving(m3_out0_mean, replace)

clear
use m3_out0_mean
twoway (bar _margin _at2 if (_at1==0)&(_at2==0)|(_at1==0)&(_at2==1)|(_at1==0)&(_at2==2), sort color(gs11) barw(.5)) ///
    (bar _margin _at2 if (_at1==0)&(_at2==3)|(_at1==0)&(_at2==4)|(_at1==0)&(_at2==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub _at2 if (_at1==0), sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel(0 "No Rel. Lev." 1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(0(.2)1) ///
	xtitle(" ") ///
	ytitle("Pr(Failure)" " ", size(large)) ///
	graphregion(color(white)) ///
	title(GDP Ratio=0, size(medlarge)) ///
	legend(order(1 2) label(1 "Without Staying Power") label(2 "With Staying Power") cols(2) region(lcolor(white)))  ///
	saving(m3_out0_gdp0, replace)

twoway (bar _margin _at2 if (_at1==.15)&(_at2==0)|(_at1==.15)&(_at2==1)|(_at1==.15)&(_at2==2), sort color(gs11) barw(.5)) ///
    (bar _margin _at2 if (_at1==.15)&(_at2==3)|(_at1==.15)&(_at2==4)|(_at1==.15)&(_at2==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub _at2 if (_at1==.15), sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel(0 "No Rel. Lev." 1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(0(.2)1) ///
	xtitle(" ") ///
	ytitle(" " " ", size(large)) ///
	graphregion(color(white)) ///
	title(GDP Ratio=Mean, size(medlarge)) ///
	legend(order(1 2) label(1 "Without Staying Power") label(2 "With Staying Power") cols(2) region(lcolor(white)))  ///
	saving(m3_out0_gdp153, replace)

twoway (bar _margin _at2 if (_at1==.27)&(_at2==0)|(_at1==.27)&(_at2==1)|(_at1==.27)&(_at2==2), sort color(gs11) barw(.5)) ///
    (bar _margin _at2 if (_at1==.27)&(_at2==3)|(_at1==.27)&(_at2==4)|(_at1==.27)&(_at2==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub _at2 if (_at1==.27), sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel(0 "No Rel. Lev." 1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(0(.2)1) ///
	xtitle(" ") ///
	ytitle(" " " ", size(large)) ///
	graphregion(color(white)) ///
	title(GDP Ratio=90 Perc., size(medlarge)) ///
	legend(order(1 2) label(1 "Without Staying Power") label(2 "With Staying Power") cols(2) region(lcolor(white)))  ///
	saving(m3_out0_gdp273, replace)
	
grc1leg m3_out0_gdp0.gph m3_out0_gdp153.gph m3_out0_gdp273.gph, saving(m3_out0_bar.gph, replace) ///
	title("") ///
		graphregion(color(white)) pos(12) cols(3) 

clear
use m3_out2_mean
twoway (bar _margin _at2 if (_at1==0)&(_at2==0)|(_at1==0)&(_at2==1)|(_at1==0)&(_at2==2), sort color(gs11) barw(.5)) ///
    (bar _margin _at2 if (_at1==0)&(_at2==3)|(_at1==0)&(_at2==4)|(_at1==0)&(_at2==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub _at2 if (_at1==0), sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel(0 "No Rel. Lev." 1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(0(.2)1) ///
	xtitle(" ") ///
	ytitle("Pr(Peace Agreement)" " ", size(large)) ///
	graphregion(color(white)) ///
	title(GDP Ratio=0, size(medlarge)) ///
    legend(off) ///
	saving(m3_out2_gdp0, replace)

twoway (bar _margin _at2 if (_at1==.15)&(_at2==0)|(_at1==.15)&(_at2==1)|(_at1==.15)&(_at2==2), sort color(gs11) barw(.5)) ///
    (bar _margin _at2 if (_at1==.15)&(_at2==3)|(_at1==.15)&(_at2==4)|(_at1==.15)&(_at2==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub _at2 if (_at1==.15), sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel(0 "No Rel. Lev." 1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(0(.2)1) ///
	xtitle(" ") ///
	ytitle(" " " ", size(large)) ///
	graphregion(color(white)) ///
	title(GDP Ratio=Mean, size(medlarge)) ///
    legend(off) ///
	saving(m3_out2_gdp153, replace)

twoway (bar _margin _at2 if (_at1==.27)&(_at2==0)|(_at1==.27)&(_at2==1)|(_at1==.27)&(_at2==2), sort color(gs11) barw(.5)) ///
    (bar _margin _at2 if (_at1==.27)&(_at2==3)|(_at1==.27)&(_at2==4)|(_at1==.27)&(_at2==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub _at2 if (_at1==.27), sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel(0 "No Rel. Lev." 1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(0(.2)1) ///
	xtitle(" ") ///
	ytitle(" " " ", size(large)) ///
	graphregion(color(white)) ///
	title(GDP Ratio=90 Perc., size(medlarge)) ///
    legend(off) ///
	saving(m3_out2_gdp273, replace)
	
graph combine m3_out2_gdp0.gph m3_out2_gdp153.gph m3_out2_gdp273.gph, saving(m3_out2_bar.gph, replace) ///
	note("Note: 95% Confidence Intervals reported.", size(medium)) ///
		graphregion(color(white))  cols(3)
		
********* Figure 4 ********
use "Complementary Replication Data.dta" 

ologit outcome3 c.gdp_ratio_thousands##i.cat_lev numofmed i.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)
margins, dydx(cat_lev) at(gdp_ratio=(0 .15 .17)) vsquish predict(outcome(2)) saving(m3_out2_dydx_cat, replace)
margins, dydx(gdp_ratio) at(cat_lev=(0(1)5)) vsquish predict(outcome(2)) saving(m3_out2_dydx_gdp, replace)
margins, dydx(cat_lev) at(gdp_ratio=(0 .15 .17)) vsquish predict(outcome(0)) saving(m3_out0_dydx_cat, replace)
margins, dydx(gdp_ratio) at(cat_lev=(0(1)5)) vsquish predict(outcome(0)) saving(m3_out0_dydx_gdp, replace)

clear
use m3_out0_dydx_cat
decode _deriv, gen(cat)
generate str1 catstr= substr(cat, 1,1)
drop cat
destring catstr, gen(cat)


twoway (bar _margin cat if (_at==1)&(cat==1)|(_at==1)&(cat==2), sort color(gs11) barw(.5)) ///
    (bar _margin cat if (_at==1)&(cat==3)|(_at==1)&(cat==4)|(_at==1)&(cat==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub cat if (_at==1), sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel(1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(-.2(.1)0.4) ///
	xtitle(" ") ///
	ytitle("Change in Pr(Failure)" " ", size(large)) ///
	graphregion(color(white)) ///
	title(GDP Ratio==0, size(medlarge)) ///
	yline(0) ///
	legend(order(1 2) label(1 "Without Staying Power") label(2 "With Staying Power") cols(2) region(lcolor(white)))  ///
	saving(m3_out0_dydx_cat_gdp0, replace)


twoway (bar _margin cat if (_at==2)&(cat==1)|(_at==2)&(cat==2), sort color(gs11) barw(.5)) ///
    (bar _margin cat if (_at==2)&(cat==3)|(_at==2)&(cat==4)|(_at==2)&(cat==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub cat if (_at==2), sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel( 1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(-.2(.1)0.4) ///
	xtitle(" ") ///
	ytitle(" " " ", size(large)) ///
	graphregion(color(white)) ///
	title(GDP Ratio==Mean, size(medlarge)) ///
	yline(0) ///
	legend(order(1 2) label(1 "Without Staying Power") label(2 "With Staying Power") cols(2) region(lcolor(white)))  ///
	saving(m3_out0_dydx_cat_gdp153, replace)

twoway (bar _margin cat if (_at==3)&(cat==1)|(_at==3)&(cat==2), sort color(gs11) barw(.5)) ///
    (bar _margin cat if (_at==3)&(cat==3)|(_at==3)&(cat==4)|(_at==3)&(cat==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub cat if (_at==3), sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel(1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(-.2(.1)0.4) ///
	xtitle(" ") ///
	ytitle(" " " ", size(large)) ///
	graphregion(color(white)) ///
	title(GDP Ratio==90 Perc., size(medlarge)) ///
		yline(0) ///
	legend(order(1 2) label(1 "Without Staying Power") label(2 "With Staying Power") cols(2) region(lcolor(white)))  ///
	saving(m3_out0_dydx_cat_gdp273, replace)
	
grc1leg m3_out0_dydx_cat_gdp0.gph m3_out0_dydx_cat_gdp153.gph m3_out0_dydx_cat_gdp273.gph, saving(m2_out0_dydx_cat.gph, replace) ///
	title("") ///
		graphregion(color(white)) pos(12)  cols(3) 		


clear
use m3_out2_dydx_cat
decode _deriv, gen(cat)
generate str1 catstr= substr(cat, 1,1)
drop cat
destring catstr, gen(cat)


twoway (bar _margin cat if (_at==1)&(cat==1)|(_at==1)&(cat==2), sort color(gs11) barw(.5)) ///
    (bar _margin cat if (_at==1)&(cat==3)|(_at==1)&(cat==4)|(_at==1)&(cat==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub cat if (_at==1), sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel(1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(-.4(.2)0.6) ///
	xtitle(" ") ///
	ytitle("Change in Pr(Peace Agreement)" " ", size(large)) ///
	graphregion(color(white)) ///
	title(GDP Ratio==0, size(medlarge)) ///
	yline(0) ///
    legend(off) ///
	saving(m3_out2_dydx_cat_gdp0, replace)


twoway (bar _margin cat if (_at==2)&(cat==1)|(_at==2)&(cat==2), sort color(gs11) barw(.5)) ///
    (bar _margin cat if (_at==2)&(cat==3)|(_at==2)&(cat==4)|(_at==2)&(cat==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub cat if (_at==2), sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel( 1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(-.4(.2)0.6) ///
	xtitle(" ") ///
	ytitle(" " " ", size(large)) ///
	graphregion(color(white)) ///
	title(GDP Ratio==Mean, size(medlarge)) ///
	yline(0) ///
    legend(off) ///
	saving(m3_out2_dydx_cat_gdp153, replace)

twoway (bar _margin cat if (_at==3)&(cat==1)|(_at==3)&(cat==2), sort color(gs11) barw(.5)) ///
    (bar _margin cat if (_at==3)&(cat==3)|(_at==3)&(cat==4)|(_at==3)&(cat==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub cat if (_at==3), sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel(1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(-.4(.2)0.6) ///
	xtitle(" ") ///
	ytitle(" " " ", size(large)) ///
	graphregion(color(white)) ///
	title(GDP Ratio==90 Perc., size(medlarge)) ///
		yline(0) ///
    legend(off) ///
	saving(m3_out2_dydx_cat_gdp273, replace)
	
graph combine m3_out2_dydx_cat_gdp0.gph m3_out2_dydx_cat_gdp153.gph m3_out2_dydx_cat_gdp273.gph, saving(m3_out2_dydx_cat.gph, replace) ///
	note("Note: Figure depicts change in predicted probability when moving from" ///
	"	   no soft leverage to each category, conditional on Relative GDP.", size(medium)) ///
		graphregion(color(white)) cols(3) 
		

********* Figure 5 ********
clear
use m3_out2_dydx_gdp
twoway (bar _margin _at2 if (_at2==0)|(_at2==1)|(_at2==2), sort color(gs11) barw(.5)) ///
    (bar _margin _at2 if (_at2==3)|(_at2==4)|(_at2==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub _at2 , sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel(0 "No Rel. Lev." 1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(-1(0.25)1.75) ///
	xtitle(" ") ///
	ytitle("Change in Pr(Peace Agt)" " ", size(large)) ///
	yline(0) ///
	graphregion(color(white)) ///
	legend(order(1 2) label(1 "Without Staying Power") label(2 "With Staying Power") cols(2) region(lcolor(white)))  ///
	saving(m3_out2_dydx_gdp, replace)

clear
use m3_out0_dydx_gdp
twoway (bar _margin _at2 if (_at2==0)|(_at2==1)|(_at2==2), sort color(gs11) barw(.5)) ///
    (bar _margin _at2 if (_at2==3)|(_at2==4)|(_at2==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub _at2 , sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel(0 "No Rel. Lev." 1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(-1(0.25)1.75) ///
	xtitle(" ") ///
	ytitle("Change in Pr(Failure)" " ", size(large)) ///
	yline(0) ///
	graphregion(color(white)) ///
	legend(order(1 2) label(1 "Without Staying Power") label(2 "With Staying Power") cols(2) region(lcolor(white)))  ///
	saving(m3_out0_dydx_gdp, replace)
	
grc1leg m3_out0_dydx_gdp.gph m3_out2_dydx_gdp.gph, saving(m3_dydx_gdp.gph, replace) ///
	note("Note: Figure depicts change in predicted probability for a one unit change in" ///
	"	   GDP ratio, conditional on soft forms of leverage.", size(medium)) ///
		graphregion(color(white)) pos(12)

		
********* Figure 6 ********
		
use "Complementary Replication Data.dta", clear		
ologit outcome3 c.gdp_ratio_thousands i.cat_lev numofmed i.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)
margins, at(coord_single=(0 1 2) numofmed=(1 3 6)) vsquish predict(outcome(2)) saving(m3_out2_coord, replace)
margins, at(coord_single=(0 1 2) numofmed=(1 3 6)) vsquish predict(outcome(0)) saving(m3_out0_coord, replace)

clear
use m3_out2_coord
drop if (_at==2)|(_at==3)|(_at==4)|(_at==7)
gen obs=_n
recode obs(3=4)(4=3)
twoway (bar _margin obs, sort color(gs10) barw(.5)) ///
    (rcap _ci_lb _ci_ub obs, sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel(1 "Single Party" 2 "Uncoord. 3 Mediators" ///
	3 "Uncoord. 6 Mediators" 	4 "Coord. 3 Mediators" ///
	5 "Coord. 6 Mediators", labsize(large) angle(45)) ///
	ylabel(0(.2)1) ///
	xtitle(" ") ///
	ytitle("Pr(Peace Agreement)" " ", size(large)) ///
	graphregion(color(white)) ///
	legend(off)  ///
	saving(m3_out2_coord, replace)

clear
use m3_out0_coord
drop if (_at==2)|(_at==3)|(_at==4)|(_at==7)
gen obs=_n
recode obs(3=4)(4=3)
twoway (bar _margin obs, sort color(gs10) barw(.5)) ///
    (rcap _ci_lb _ci_ub obs, sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel(1 "Single Party" 2 "Uncoord. 3 Mediators" ///
	3 "Uncoord. 6 Mediators" 	4 "Coord. 3 Mediators" ///
	5 "Coord. 6 Mediators", labsize(large) angle(45)) ///
	ylabel(0(.2)1) ///
	xtitle(" ") ///
	ytitle("Pr(Failure)" " ", size(large)) ///
	graphregion(color(white)) ///
	legend(off)  ///
	saving(m3_out0_coord, replace)

	
graph combine m3_out0_coord.gph m3_out2_coord.gph, saving(m3_coord_bar.gph, replace) ///
	note("Note: 95% Confidence Intervals reported.", size(medium)) ///
		graphregion(color(white)) 
		
***************************************
************* Footnotes ***************
***************************************
use "Complementary Replication Data.dta", clear		

***** Potential Multicolinearity *****
reg outcome3 c.gdp_ratio_thousands i.cat_lev numofmed i.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)
vif	
reg outcome3 c.gdp_ratio_thousands##i.cat_lev numofmed i.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)
vif

reg outcome3 c.gdp_ratio_thousands i.cat_lev numofmed i.coord_single fullconf_duration_10years intensity_level terr_conflict currentpk polity2_v5, cl(conflictid)
vif 
reg outcome3 c.gdp_ratio_thousands##i.cat_lev numofmed i.coord_single fullconf_duration_10years intensity_level terr_conflict currentpk polity2_v5, cl(conflictid)
vif

pwcorr gdp_ratio_thousands cat_lev numofmed coord_single fullconf_duration_10years currentpk intensity_level terr_conflict polity2_v5, sig

***** Current v Promise of PKing *****
tab stay_any currentpk, row chi2




***************************************
************* Appendix ****************
***************************************


***** Alternative Operationalizations *****
use "Complementary Replication Data.dta", clear		

** Table A1. Log of GDP max instead of GDP ratio 
*gen lnmax_gdp_med = ln(max_gdp_med+1)
ologit outcome3 c.lnmax_gdp_med i.cat_lev numofmed i.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)
estimates store m1
ologit outcome3 c.lnmax_gdp_med##i.cat_lev numofmed i.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)
estimates store m2

** Table A1. CINC ratio instead of GDP ratio 
ologit outcome3 c.cinc_ratio i.cat_lev numofmed i.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)
estimates store m3
ologit outcome3 c.cinc_ratio##i.cat_lev numofmed i.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)
estimates store m4

** Table A1. Log of Num of Med 
ologit outcome3 c.gdp_ratio_thousands i.cat_lev log_med i.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)
estimates store m5
ologit outcome3 c.gdp_ratio_thousands##i.cat_lev log_med i.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)
estimates store m6

ologit outcome3 c.gdp_ratio_thousands i.cat_lev log_med ib1.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)
ologit outcome3 c.gdp_ratio_thousands##i.cat_lev log_med ib1.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)

esttab m1 m2 m3 m4 m5 m6 using appendix_table1.tex, b(3) compress ///
		se label replace cells(b(fmt(3) star)  ///
		se(par)) starlevels(* .05)  ///
		varwidth(30) collabels(, none) noomitted
		
************ Figure A1 ***************
ologit outcome3 c.lnmax_gdp_med##i.cat_lev numofmed i.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)
margins, dydx(lnmax_gdp_med) at(cat_lev=(0 1 2 3 4 5)) vsquish predict(outcome(2)) saving(m2_out2_dydx_gdp, replace)
margins, dydx(lnmax_gdp_med) at(cat_lev=(0 1 2 3 4 5)) vsquish predict(outcome(0)) saving(m2_out0_dydx_gdp, replace)

clear
use m2_out2_dydx_gdp
twoway (bar _margin _at2 if (_at2==0)|(_at2==1)|(_at2==2), sort color(gs11) barw(.5)) ///
    (bar _margin _at2 if (_at2==3)|(_at2==4)|(_at2==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub _at2 , sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel(0 "No Rel. Lev." 1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(-0.05(0.02)0.075) ///
	xtitle(" ") ///
	ytitle("Change in Pr(Peace Agt)" " ", size(large)) ///
	yline(0) ///
	graphregion(color(white)) ///
	legend(order(1 2) label(1 "Without Staying Power") label(2 "With Staying Power") cols(2) region(lcolor(white)))  ///
	saving(m2_out2_dydx_gdp, replace)

clear
use m2_out0_dydx_gdp
twoway (bar _margin _at2 if (_at2==0)|(_at2==1)|(_at2==2), sort color(gs11) barw(.5)) ///
    (bar _margin _at2 if (_at2==3)|(_at2==4)|(_at2==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub _at2 , sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel(0 "No Rel. Lev." 1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(-0.05(0.02)0.075) ///
	xtitle(" ") ///
	ytitle("Change in Pr(Failure)" " ", size(large)) ///
	yline(0) ///
	graphregion(color(white)) ///
	legend(order(1 2) label(1 "Without Staying Power") label(2 "With Staying Power") cols(2) region(lcolor(white)))  ///
	saving(m2_out0_dydx_gdp, replace)
	
grc1leg m2_out0_dydx_gdp.gph m2_out2_dydx_gdp.gph, saving(m2_dydx_gdp.gph, replace) ///
	note("Note: Figure depicts change in predicted probability for a one unit change in" ///
	"	   Ln(GDP), conditional on soft forms of leverage.", size(medium)) ///
		graphregion(color(white)) pos(12)

************ Figure A2 ***************
use "Complementary Replication Data.dta", clear		

ologit outcome3 c.cinc_ratio##i.cat_lev numofmed i.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)
margins, dydx(cinc_ratio) at(cat_lev=(0 1 2 3 4 5)) vsquish predict(outcome(2)) saving(m4_out2_dydx_gdp, replace)
margins, dydx(cinc_ratio) at(cat_lev=(0 1 2 3 4 5)) vsquish predict(outcome(0)) saving(m4_out0_dydx_gdp, replace)

clear
use m4_out2_dydx_gdp
twoway (bar _margin _at2 if (_at2==0)|(_at2==1)|(_at2==2), sort color(gs11) barw(.5)) ///
    (bar _margin _at2 if (_at2==3)|(_at2==4)|(_at2==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub _at2 , sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel(0 "No Rel. Lev." 1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(-0.003(0.002)0.007) ///
	xtitle(" ") ///
	ytitle("Change in Pr(Peace Agt)" " ", size(large)) ///
	yline(0) ///
	graphregion(color(white)) ///
	legend(order(1 2) label(1 "Without Staying Power") label(2 "With Staying Power") cols(2) region(lcolor(white)))  ///
	saving(m4_out2_dydx_gdp, replace)

clear
use m4_out0_dydx_gdp
twoway (bar _margin _at2 if (_at2==0)|(_at2==1)|(_at2==2), sort color(gs11) barw(.5)) ///
    (bar _margin _at2 if (_at2==3)|(_at2==4)|(_at2==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub _at2 , sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel(0 "No Rel. Lev." 1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(-0.003(0.002)0.007) ///
	xtitle(" ") ///
	ytitle("Change in Pr(Failure)" " ", size(large)) ///
	yline(0) ///
	graphregion(color(white)) ///
	legend(order(1 2) label(1 "Without Staying Power") label(2 "With Staying Power") cols(2) region(lcolor(white)))  ///
	saving(m4_out0_dydx_gdp, replace)
	
grc1leg m4_out0_dydx_gdp.gph m4_out2_dydx_gdp.gph, saving(m4_dydx_gdp.gph, replace) ///
	note("Note: Figure depicts change in predicted probability for a one unit change in" ///
	"	   CINC ratio, conditional on soft forms of leverage.", size(medium)) ///
		graphregion(color(white)) pos(12)

			
		
***** Alternative Operationalization of DV *****
** Table A2. 4 Category Outcome as DV
use "Complementary Replication Data.dta", clear		

ologit outcome4 c.gdp_ratio_thousands i.cat_lev numofmed i.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)
estimates store m5
ologit outcome4 c.gdp_ratio_thousands##i.cat_lev numofmed i.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)
estimates store m6

esttab m5 m6 using appendix_table3.tex, b(3) compress ///
		se label replace cells(b(fmt(3) star)  ///
		se(par)) starlevels(* .05)  ///
		varwidth(30) collabels(, none) noomitted


***** Controlling for Reputation *****
** Table A3. Including Reputation 

ologit outcome3 c.gdp_ratio_thousands i.cat_lev rep_prop numofmed i.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)
estimates store rep
ologit outcome3 c.gdp_ratio_thousands##i.cat_lev rep_prop numofmed i.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)
estimates store rep2

esttab rep rep2 using appendix_table3.tex, b(3) compress ///
		se label replace cells(b(fmt(3) star)  ///
		se(par)) starlevels(* .05)  ///
		varwidth(30) collabels(, none) noomitted
		
ologit outcome3 c.gdp_ratio_thousands i.cat_lev rep_any numofmed i.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)
ologit outcome3 c.gdp_ratio_thousands##i.cat_lev rep_any numofmed i.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)

** Figure A3. Including Reputation 

ologit outcome3 c.gdp_ratio_thousands##i.cat_lev rep_prop numofmed i.coord_single fullconf_duration_10years intensity_level terr_conflict polity2_v5, cl(conflictid)
margins, dydx(cat_lev) at(gdp_ratio=(0 .15 .17)) vsquish predict(outcome(2)) saving(repmed, replace)
margins, dydx(cat_lev) at(gdp_ratio=(0 .15 .17)) vsquish predict(outcome(0)) saving(repmed_fail, replace)

clear
use repmed_fail
decode _deriv, gen(cat)
generate str1 catstr= substr(cat, 1,1)
drop cat
destring catstr, gen(cat)


twoway (bar _margin cat if (_at==1)&(cat==1)|(_at==1)&(cat==2), sort color(gs11) barw(.5)) ///
    (bar _margin cat if (_at==1)&(cat==3)|(_at==1)&(cat==4)|(_at==1)&(cat==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub cat if (_at==1), sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel(1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(-.2(.1)0.4) ///
	xtitle(" ") ///
	ytitle("Change in Pr(Failure)" " ", size(large)) ///
	graphregion(color(white)) ///
	title(GDP Ratio==0, size(medlarge)) ///
	yline(0) ///
	legend(order(1 2) label(1 "Without Staying Power") label(2 "With Staying Power") cols(2) region(lcolor(white)))  ///
	saving(m3_out0_dydx_cat_gdp0, replace)


twoway (bar _margin cat if (_at==2)&(cat==1)|(_at==2)&(cat==2), sort color(gs11) barw(.5)) ///
    (bar _margin cat if (_at==2)&(cat==3)|(_at==2)&(cat==4)|(_at==2)&(cat==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub cat if (_at==2), sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel( 1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(-.2(.1)0.4) ///
	xtitle(" ") ///
	ytitle(" " " ", size(large)) ///
	graphregion(color(white)) ///
	title(GDP Ratio==Mean, size(medlarge)) ///
	yline(0) ///
	legend(order(1 2) label(1 "Without Staying Power") label(2 "With Staying Power") cols(2) region(lcolor(white)))  ///
	saving(m3_out0_dydx_cat_gdp153, replace)

twoway (bar _margin cat if (_at==3)&(cat==1)|(_at==3)&(cat==2), sort color(gs11) barw(.5)) ///
    (bar _margin cat if (_at==3)&(cat==3)|(_at==3)&(cat==4)|(_at==3)&(cat==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub cat if (_at==3), sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel(1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(-.2(.1)0.4) ///
	xtitle(" ") ///
	ytitle(" " " ", size(large)) ///
	graphregion(color(white)) ///
	title(GDP Ratio==90 Perc., size(medlarge)) ///
		yline(0) ///
	legend(order(1 2) label(1 "Without Staying Power") label(2 "With Staying Power") cols(2) region(lcolor(white)))  ///
	saving(m3_out0_dydx_cat_gdp273, replace)
	
grc1leg m3_out0_dydx_cat_gdp0.gph m3_out0_dydx_cat_gdp153.gph m3_out0_dydx_cat_gdp273.gph, saving(m2_out0_dydx_cat.gph, replace) ///
	title("") ///
		graphregion(color(white)) pos(12)  cols(3) 		


clear
use repmed
decode _deriv, gen(cat)
generate str1 catstr= substr(cat, 1,1)
drop cat
destring catstr, gen(cat)


twoway (bar _margin cat if (_at==1)&(cat==1)|(_at==1)&(cat==2), sort color(gs11) barw(.5)) ///
    (bar _margin cat if (_at==1)&(cat==3)|(_at==1)&(cat==4)|(_at==1)&(cat==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub cat if (_at==1), sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel(1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(-.4(.2)0.6) ///
	xtitle(" ") ///
	ytitle("Change in Pr(Peace Agreement)" " ", size(large)) ///
	graphregion(color(white)) ///
	title(GDP Ratio==0, size(medlarge)) ///
	yline(0) ///
    legend(off) ///
	saving(m3_out2_dydx_cat_gdp0, replace)


twoway (bar _margin cat if (_at==2)&(cat==1)|(_at==2)&(cat==2), sort color(gs11) barw(.5)) ///
    (bar _margin cat if (_at==2)&(cat==3)|(_at==2)&(cat==4)|(_at==2)&(cat==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub cat if (_at==2), sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel( 1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(-.4(.2)0.6) ///
	xtitle(" ") ///
	ytitle(" " " ", size(large)) ///
	graphregion(color(white)) ///
	title(GDP Ratio==Mean, size(medlarge)) ///
	yline(0) ///
    legend(off) ///
	saving(m3_out2_dydx_cat_gdp153, replace)

twoway (bar _margin cat if (_at==3)&(cat==1)|(_at==3)&(cat==2), sort color(gs11) barw(.5)) ///
    (bar _margin cat if (_at==3)&(cat==3)|(_at==3)&(cat==4)|(_at==3)&(cat==5), sort color(gs7) barw(.5)) ///
    (rcap _ci_lb _ci_ub cat if (_at==3), sort pstyle(ci) color(black) lw(medthick)), ///
    xlabel(1 "Rel. 1 Side" ///
	2 "Rel. Both Sides" 	3 "No Rel. Lev." ///
	4 "Rel. 1 Side" 5 "Rel. Both Sides", labsize(large) angle(45)) ///
	ylabel(-.4(.2)0.6) ///
	xtitle(" ") ///
	ytitle(" " " ", size(large)) ///
	graphregion(color(white)) ///
	title(GDP Ratio==90 Perc., size(medlarge)) ///
		yline(0) ///
    legend(off) ///
	saving(m3_out2_dydx_cat_gdp273, replace)
	
graph combine m3_out2_dydx_cat_gdp0.gph m3_out2_dydx_cat_gdp153.gph m3_out2_dydx_cat_gdp273.gph, saving(m3_out2_dydx_cat.gph, replace) ///
	note("Note: Figure depicts change in predicted probability when moving from" ///
	"	   no soft leverage to each category, conditional on Relative GDP.", size(medium)) ///
		graphregion(color(white)) cols(3) 
		

***** Model Selection *****
use "Complementary Replication Data.dta", clear		

** Table A4. Gologit2 
*gologit doesn't like factor variables; creating dummies
tab cat_lev, g(cat_lev_dummies)
tab coord_single, g(coord_single_dummies)

gologit2 outcome3 gdp_ratio_thousands cat_lev_dummies2 cat_lev_dummies3 cat_lev_dummies4 cat_lev_dummies5 cat_lev_dummies6 numofmed coord_single_dummies2 coord_single_dummies3 fullconf_duration_10years intensity_level terr_conflict polity2_v5, autofit cl(conflictid)
estimates store m9

*gologit doesn't like interactions; creating them here
gen gdp_cat_1= gdp_ratio_thousands * cat_lev_dummies2
gen gdp_cat_2= gdp_ratio_thousands * cat_lev_dummies3
gen gdp_cat_3= gdp_ratio_thousands * cat_lev_dummies4
gen gdp_cat_4= gdp_ratio_thousands * cat_lev_dummies5
gen gdp_cat_5= gdp_ratio_thousands * cat_lev_dummies6

gologit2 outcome3 gdp_ratio_thousands cat_lev_dummies2 cat_lev_dummies3 cat_lev_dummies4 cat_lev_dummies5 cat_lev_dummies6 gdp_cat_1 gdp_cat_2 gdp_cat_3 gdp_cat_4 gdp_cat_5 numofmed coord_single_dummies2 coord_single_dummies3 fullconf_duration_10years intensity_level terr_conflict polity2_v5, autofit cl(conflictid)
estimates store m10

esttab m9 m10 using appendix_table4.tex, b(3) compress ///
		se label replace cells(b(fmt(3) star)  ///
		se(par)) starlevels(* .05)  ///
		varwidth(30) collabels(, none) noomitted
		
	
** Table A5. Selection Models
use "Complementary ConflictMonth Replication Data.dta", clear	
heckoprobit outcome3 c.gdp_ratio_thousands i.cat_lev numofmed i.coord_single conf_dur_selection_10yr intensity_level incompatibility_ac_terr polity2_v5 , select(newmed=parity lag_mediation incompatibility_ac_terr conf_dur_selection_10yr) cl(conflict_id)
estimates store sel1
heckoprobit outcome3 c.gdp_ratio_thousands##i.cat_lev numofmed i.coord_single conf_dur_selection_10yr intensity_level incompatibility_ac_terr polity2_v5 , select(newmed=parity lag_mediation incompatibility_ac_terr conf_dur_selection_10yr) cl(conflict_id)
estimates store sel2


esttab sel1 sel2 using appendix_table5.tex, b(3) compress ///
		se label replace cells(b(fmt(3) star)  ///
		se(par)) starlevels(* .05)  ///
		varwidth(30) collabels(, none) noomitted
